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Abstract: Electromagnetic inverse scattering problems for point-like scatterers are considered in 
this presentation. A scatterer can be treated as point-like when its size is much smaller than one 
wavelength. Although the inverse scattering problem is a nonlinear problem, it can be solved 
without iteratively evaluating the corresponding forward problem, which means that the problem 
can be solved very fast. In this presentation, the multiple signal classification (MUSIC) method 
and the two-step least squares method are provided, with the former being used to locate point- 
like scatterers and the latter being used to retrieve the scattering strengths of point-like scatterers. 
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Description of the problem 


= Determine the locations and polarization tensors (scattering 
coefficient tensors) of a collection of small scatterers from 
scattered data (multistatic response matrix) 
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Forward scattering problem (General Case) 


= Homogeneous background with & % 
= Distribution of permittivities <.(7) 
= Incidence field p" : 


VxVxE -kE =0 (a) 


wave number: k=a,/é 4 
= Total field g” : 
we lk ga — tot = —tot 
VxE =i@H VxH =-iwé,é-E 


— tot 


VxVxE -ke -E =0 (b) 
= Scattered field 


— sca — tot —inc 


E =E - 





Forward scattering problem (General Case) 


= Scattered field is due to induced current: 


— tot 


VxVxE OKRE” tire Be 


J = E-E” , Where é depends on E, 
= From Green’s function, we have in the domain of interest D 


—sca — — — — —d — — 
E (r)=|[ Grr) J (rddr' 
= Lippmann-Schwinger equation: for reD 
E (r)=E (r)+ f [Grr (ddr 


= The relationship between incidence field and measured scattered 
field: o _ 
E” (Free) _ K.E" (Taom) 


= To reconstruct E(r) from kK isan ill-posed, nonlinear, and high- 
dimensional problem. 
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Forward scattering problem 


= Antennas: 
> Total number of sets: N 
> Each set consists of three antennas in the x, y, z directions. 


: / / / 
> Locations: T4;V'o,---,Ty 


= Scatterers: 
> Total number: M 
> Much smaller than wavelength à: ka «1 
> Locations: 1, ť2,..., rM 


> Polarization tensors: E} 


Il(r;) = £; - Ett (r;) 





Forward scattering problem 


= The leading term in scattered field is on the order of (ka) 


"= Lippmann-Schwinger equation is reduced to Foldy-Lax equation 


= Multistatic response (MSR) matrix 
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RG, j) = twpoG(r4, ry) 


= diag £1, €9, es „ÊM. 
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Inverse Scattering: Overview 


Consider two questions: 


= Is it possible to solve this inverse problem in an approach 
that does not iteratively evaluate the corresponding forward 
scattering problem’? 
> Positioning 
> Retrieval of polarization tensors 
= How to optimally implement the algorithm? 








= For small scatterers: The mapping from currents in finite 
discrete scatterers to scattered fields is one-to-one. 


(Scalar wave: A. Kirsch: Inverse Problems, 2002, Theorem 2.1 ) 


> Green's function evaluated at a point cannot be expressed as a 
linear combination of Green's functions evaluated at other discrete 
points. — 

G(r, ro) aQ £ 4 G(r, r;) ai 
= Positions can be obtained by Multiple-Signal-Classification (MUSIC) 
> Induced currents can be uniquely obtained from scattered fields 
SM G(r, r;)a; =O iff aj=O,2=—1,2,...,M 
= Polarization tensors can be obtained by least-squares method 


= In the whole inverse process, the corresponding forward 
scattering problem Is not iteratively evaluated. 
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Inverse Scattering: MUSIC 


=" MUSIC applied to detect objects: 


> In acoustics: 
= A. J. Devaney, 2000 
= S. Hou, /nverse Problems, 2006, 2007 
> A special case of factorization method 
= A. Kirsch: Inverse Problems, 2002 
= M. Cheney: Inverse Problems, 2001 
> In 3D electromagnetics: 
= H. Ammari: SIAM, 2007 
= Y. Zhong, IEEE Trans. AP, 2007 
> Phase-free 
= E. A. Marengo, JOSA, A, 2007 
= X. Chen, JOSA, A, 2008 





Standard MUSIC algorithm 


= The vector of scattered fields at receivers is in the space spanned 
by the background Green's function vectors associated with the x, 
y, and z components of electric dipoles evaluated the position of 
each scatterer 


ES € So = span {Gr(r;), Gy(r;), Ge(rj); f = 1,2,...,M} 


= Singular value decomposition: 


K . O38N — S, c C3N 
C3N = SO Sn 
Sr = span{up, cp > O} 


Sn = Span{up, op = O} 


= S; and So are identical a401 (rm) = 0 op = O 
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Standard MUSIC algorithm 


" Define MUSIC pseudo-spectrum 
P(r) — — 
Dool E) 
f(r) is any linear combination of 
G;(r), G,(r), G(r) 


The pseudo-spectrum becomes infinite at the position of 
each scatterer 


= In electromagnetic inverse scattering problem, standard 
MUSIC is a direct generalization of MUSIC in acoustics. 

= It assumes that there are three freedoms in the induced 
dipole at the location of each scatterer. 





Limitations of standard MUSIC 





= The standard MUSIC is not capable of dealing with 
degenerate scatterers, 


> needle- or disk-shaped scatterers 
> scatterers with special anisotropic material (e.g. E€,= £o). 


Mathematically, degenerate scatterers: a rank-deficient 
polarization tensor. 


Isotropic 


= In EM scattering scenario, the Sphere 
most fundamental induced 


source is dipole, unlike the one Disk \ 
es 





(monopole) in acoustic imaging. 


> The choose of the direction of test Sf 


dipole determines imaging 2 
resolution. 
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A New MUSIC Algorithm 


= Aiming at solving these two problems in a unified scheme: 
1) To detect degenerate scatterers 
2) To enhance resolution 


= Natural idea: The problem of detecting degenerate scatterers 
is solved by choosing the direction of test dipole that 
minimizes the projection of the test Green's function vector 
onto the noise space. 


= However, it cannot be used to enhance the resolution due to 
the intrinsic instability of the noise subspace. Need to avoid 
the involvement of the noise subspace in a MUSIC algorithm. 

= Since the range of the MSR are spanned by those singular 
vectors associated to leading singular values, we thus turn to 
the range, i.e. to find a test dipole direction that produces a 
Green's function vector robustly belonging to the range. 


Details can be found in X. Chen et.al., Inverse Problems, V.25, 015008 (12pp), 2009 
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A New MUSIC Algorithm 


= We try to construct a stable test direction by using 

singular vectors in the range 
S4] Aŭ; = G(r) -a Jall = 1 and a € C? 

= Meanwhile, the singular vector can also be represented 
by the Green's function vector generated at the positions 
of the scatterers 

a; = OM, Gry) -IM,  §=1,2,...,L 

where J” is the current source induced inside the j" 
scatterer in the i eigen-state . 

= Use the independency of Green's function vectors. 


e.g. when the test point is at the position of the first 
scatterer 
DiJ Pasa DJP =0, j=2,3,...M 
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A New MUSIC Algorithm 


= Only one dipole at the test point is left, and dipoles at all 
other scatterers are zero. Thus, L is equal to one plus 
the total independent dipoles induced at all other dipoles. 
> e.g. M isotropic spherical scatterers, L = 3(M-1)+7, whereas the 
rank of MSR matrix is 3M. 
= In practice, L is determined from the dropping trend of 
the singular values of the MSR matrix. 


a EAUS; the direction of the test dipole could be determined 


y _, |ut-G(r)-al? 
amax = arg max, atti Gu) al’ iC (r)al 


|G(r)-al? 
= And anew pseudo-spectrum could be defined as 
b(r) = —~—+-.—_, 
( ) | oie fC @)-amax|? 
|G(r)-amax|? 


= The a obtained is generally complex. If we constrain a to 
be real (i.e., physical spatial direction), we will need a 
larger L, thus losing some stability. 





of Singapore 





Inverse Scattering: Overview 


= Retrieval of scattering strength 
> Iterative numerical approach: 
= A. J. Devaney, J. Acoust. Soc. Am., 2005 
> Non-iterative analytical approach 
= E. A. Marengo, J. Acoust. Soc. Am., 2006 
= X. Chen, J. Acoust. Soc. Am., 2007, 
= Y. Zhong, IEEE, Tans. AP, 2007 
> Phase-free 
= X. Chen, , JOSA, A, 2008 
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Inverse Scattering: Least Squares 


= Retrieve polarization tensors from MSR matrix using a 
non-iterative analytical method. 


K=R-A-(I3y—-®-A)7!- 


» The induced current is represente 
J = — A. Gomu- P. A)~ 


= 
L = = 


its least squares solution 
J=R-K 
> New equation (Foldy-Lax equation): 
J=A-(F+-J) 
J=: (T+. J)) i = 1,2,... 
> Least squares solution 


R=, [(F+3-7)) 
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= Three scatterers are placed around the origin: 
r, = (0.084X, 0.196, 0.0842) €y = €9 = 2€ 
ro = (—0.168, —0.056\, —0.112)) Eq = diag [é9, 3€0, 9€p) 


r3 = (—0.196A, —0.084\, 0.140)) (1), 6,0) = (4/4, 7/3, 37/8) 


lmin = 0.255A y = x + 0.112A 
= 16 antenna units, half of which are aligned along the y 
axis while the other half aligned along the z axis in the 
x = —13, plane. The two linear arrays are centered at 


(—13\,—-9\,11\) with 5A separation distance between 
neighboring units. 





23 


Simulations — Without Noise 












o 0.1 
x (A) 
(c) 


Singular values and pseudo-spectrum obtained by the 
standard MUSIC algorithm in noise free case. 
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Simulations — Without Noise 
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Pseudo-spectrum obtained by the proposed MUSIC algorithm 
in noise free case. 





20 


Simulations — With Noise 
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“6 
Singular values and pseudo-spectrum obtained by the 
standard MUSIC algorithm in noise-contaminated case (3%). 
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Simulations — With Noise 
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=a Pseudo-spectrum 


obtained by the proposed 
MUSIC algorithm in 
noise-contaminated case 
(3%). 





1500 


1000 


2/ 


"Inverse scattering problem of determining the locations and the 
polarizations of small scatterers 1s solved by a non-iterative 
analytical approach. 


"MUSIC method 1s used to locate objects. 
«Least squares method is used to obtain the polarization tensors. 


"Manipulating on the range, instead of noise space, of the MSR 
matrix, we can solving these two problems in a unified scheme: 
|) To detect degenerate scatterers; 
2) To enhance resolution. 
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